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ABSTRACT 

Numerous spectroscopic and photometric observational studies have provided strong 
evidence for the widespread presence of multiple stellar populations in globular clus- 
ters. In this paper we study the long-term dynamical evolution of multiple-population 
i clusters, focusing on the evolution of the spatial distributions of the first- (FG) and 

^ second-generation (SG) stars. In previous studies we have suggested that SG stars 

formed from the ejecta of FG AGB stars are expected initially to be concentrated in 
the cluster inner regions. Here, by means of TV-body simulations, we explore the time 
^ C| scales and the dynamics of the spatial mixing of the FG and the SG populations and 

Q H their dependence on the SG initial concentration. 

Our simulations show that, as the evolution proceeds, the radial profile of the 
SG/FG number ratio, Nsg/Nfg, is characterized by three regions: 1) a flat inner part; 
2) a declining part in which FG stars are increasingly dominant; and 3) an outer region 
^ where the Nsg/^fg profile flattens again (the Nsg/^fg profile may rise slightly 

i— i again in the outermost cluster regions). Until mixing is complete and the Nsg/^fg 

profile is flat over the entire cluster, the radial variation of Nsg/^fg implies that 
r ^ the fraction of SG stars determined by observations covering a limited range of radial 

distances is not, in general, equal to the SG global fraction, (NsG/^FG)giob' The 
distance at which Nsg/^fg equals (NsGf^FG)giob is approximately between 1 and 
\^ 2 cluster half-mass radii. The time scale for complete mixing depends on the SG 

initial concentration, but in all cases complete mixing is expected only for clusters in 
advanced evolutionary phases, having lost at least 60-70 percent of their mass due 
to two-body relaxation (in addition to the early FG loss due to the cluster expansion 
triggered by SNII ejecta and gas expulsion). The results of our simulations suggest that 
in many Galactic globular clusters the SG should still be more spatially concentrated 
than the FG. 

> 

• i-H Key words: globular clustersigeneral, stars: chemically peculiar, methods:N-body 

^ simulations 



1 INTRODUCTION 

An increasing number of spectroscopic and photometric 
observational studies have provided strong evidence for a 
widespread presence of multiple stellar populations in glob- 
ular clusters. The observed star-to- star variations of light 
elements, such as Na, O, Al, and Mg, indicates that a signif- 
icant fraction (50-80%) of globular cluster stars must have 
formed out of matter processed through a high-temperature 
CNO cycle in a first generation of stars (hereafter FG; see 
e.g. Carretta et al. 2009a, 2009b and references therein). 
Photometric studies, by revealing the presence of multiple 
main sequences, subgiant, and red-giant branches in numer- 



ous clusters, have buttressed the spectroscopic evidence and 
added important new elements to the observational frame- 
work of multiple population in globular clusters. In par- 
ticular, photometric measurements have shown that some 
clusters host a population of very He-rich stars among the 
second-generation (hereafter SG) population of some clus- 
ters (see e.g. Piotto et al. 2007; see also D Antona et al. 2002, 
DAntona & Caloi 2004, D Antona & Caloi 2008). Recently, 
the first direct measurement of He abundances of two stars 
in NGC 2808 has confirmed the strong He-enhancement sug- 
gested by those photometric studies (Pasquini et al. 2011; 
see also Dupree et al. 2011 for evidence of He enhancement 
in oj Cen stars). 
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Models for the source of polluted gas from which SG 
stars formed include AGB stars (see e.g. Cottrell & Da Costa 
1981, Ventura et al. 2001), rapidly rotating massive stars 
(Decressin et al. 2007), and massive binary stars (De Mink 
et al. 2009). A few studies have also explored some aspects of 
cluster formation and evolution for some of these models (see 
e.g. D'Ercole et al. 2008, 2010, 2012, Bekki 2011, Decressin 
et al. 2008, 2010; see Renzini 2008 and Gratton et al. 2012 
and references therein for a review). 

In D'Ercole et al. (2008), we focused our attention on 
the AGB model and explored the formation and dynamical 
evolution of multiple populations by means of hydrodynam- 
ical and N-body simulations. Our simulations show that the 
AGB eject a form a cooling flow and rapidly collect in the 
innermost regions of the cluster, forming a concentrated SG 
stellar subsystem (see also Bekki 2011). In order to form the 
number of SG stars observed today, the FG cluster must 
have been considerably more massive than it is now. The N- 
body simulations presented in D'Ercole et al. (2008) show 
that the early expansion triggered by the loss of mass in the 
form of SNII ejecta leads to a strong preferential loss of FG 
stars resulting in a cluster in which the number of SG stars is 
similar to (or even larger than) that of FG stars as observed 
in several Galactic globular clusters (see e.g. Carretta et al. 
2009a, 2009b). 

In subsequent papers we expanded the initial models 
presented in D'Ercole et al. (2008) to further explore the ori- 
gin of the observed abundance patterns (D'Ercole et al. 2010, 
2012), the connection between multiple population globular 
clusters and the Galactic stellar halo (Vesperini et al. 2010) 
and the implications for the disruption of FG and SG binary 
stars (Vesperini et al. 2011). 

In this paper we focus our attention on the long- 
term dynamical evolution of the multiple-population clus- 
ter, starting immediately after the early evolutionary stages 
during which a large fraction of the FG population is lost. 
We emphasize that, throughout this paper, when we discuss 
cluster mass loss we are referring only to mass lost by the 
cluster during its long-term relaxation-driven evolution, not 
the early phase responsible for the loss of most of the initial 
FG population. 

After the early loss of FG stars (e.g. at « 1 — 2 Gyr 
for the simulations presented in D'Ercole et al. 2008), a 
multiple-population cluster will start its long-term evolution 
driven by two-body relaxation with a similar number of SG 
and FG stars but with the SG population still concentrated 
in the cluster inner regions. A number of observational stud- 
ies (Bellini et al. 2009, Carretta et al. 2010a, Lardo et al. 
2011, Kravtsov et al., 2010, 2011, Nataf et al. 2011, Johnson 
& Pilachowski 2012, Milone et al. 2012) have found that in 
several clusters SG stars are indeed preferentially located in 
the inner regions and retain some memory of the initial seg- 
regation of the SG population predicted by the models we 
presented in D'Ercole et al. (2008). 

Understanding the dynamics of the spatial mixing, and 
specifically the extent to which memory of the initial SG 
segregation predicted by the formation models is retained 
after one Hubble time of relaxation-driven evolution, is an 
essential step to properly interpreting observational data, as 
well as testing the key elements of theoretical scenarios of 
cluster formation and evolution. 

In this paper, by means of N-body simulations, we ex- 



plore the structural evolution of multiple population clus- 
ters. We focus our attention on the spatial mixing of the SG 
and the FG populations, and the evolution of the relative 
spatial distribution of the SG and the FG stars. The struc- 
ture of the paper is as follows: in Section 2 we describe the 
initial conditions of our N-body simulations; in Section 3 we 
present our results; and in Section 4 we discuss our results 
and summarize our main conclusions. 



2 METHOD AND INITIAL CONDITIONS 

The study presented in this paper is based on N-body simu- 
lations run with the starlab package (Portegies Zwart et al. 
2001) and accelerated by GRAPE-6 special purpose hard- 
ware (Makino et al. 2003) and GPU (Gaburov et al. 2009). 
We explore the evolution of four systems with different de- 
grees of initial concentration of the SG population. In all 
cases the initial conditions are set by simply combining two 
single-mass King (1966) models: the FG cluster is mod- 
eled as a King model with dimensionless central potential 
Wo = 7; the SG system is also modeled as a King model 
with Wo — 7, but it is initially entirely contained within the 
inner regions of the FG system. Our models have an initial 
ratio of the FG to the SG half-mass radius, Rh,fg / Rh,sg 
equal to 2.5, 5, 10, and 25 (Rh,FG and Rh,SG are the 3D half- 
mass radii). Hereafter we refer to these simulations as r2p5, 
r5, rlO, and r25. The models presented do not include pri- 
mordial binaries or stellar evolution. The combined system 
is initially scaled to a virial ratio T/V — 0.5 and is allowed 
to dynamically settle into an equilibrium state during the 
first few dynamical times of each simulation. 

All simulations start with a total number of particles 
N = 10, 000; initially the numbers of FG and SG particles 
are equal. We assume that the cluster is tidally truncated 
and therefore that the King truncation radius is equal to 
the cluster Jacobi radius, Rj. Particles are removed from 
the simulation after reaching a radius equal to twice the Ja- 
cobi radius. As discussed in Section [I] the initial conditions 
adopted for our simulations refer to the possible structural 
properties of a multiple-population cluster after the early 
loss of FG stars. 

To explore the role of the relaxation-driven evapora- 
tion of stars on the cluster structural and mixing properties, 
we also include a simulation with the same initial condi- 
tions as those of the r25 system, but without a tidal field or 
tidal truncation, and in which stars were removed only after 
reaching a distance of about 1000 times the cluster half-mass 
radius. Hereafter we refer to this simulation as r25isol. 

All simulations were run until the total mass of the 
system was 20 percent or less of its initial value. For the 
r25isol system, mass loss proceeds at a much slower rate 
and the simulation was stopped after the system had lost 20 
percent of its initial mass. In order to explore the dependence 
of our results on the number of particles, two additional 
simulations with N = 30,000 and N = 60,000 particles 
were also run, with initial structural properties identical to 
those of the rlO system. 
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Figure 1. Time evolution of the 1%, 5%, 10%, 25%, 50%, 75%, 90% lagrangian radii for SG stars (blue lines), FG stars (red lines) and 
all stars (black lines). The five rows of figures refer (from top to bottom) to the r25, r25isol, rlO, r5, and r2p5 simulations. In each panel, 
time is expressed relative to the initial half-mass relaxation time of the entire system in the bottom horizontal axis, and to the initial 
half-mass relaxation time of the SG subsystem in the upper horizontal axis. In order to show more clearly the details of the lagrangian 
radius evolution, for each simulation (except for r2p5) the left panels show only the evolution until ~ 25£ r ^(0) and the right panels the 
rest of the evolution until the two populations are mixed. 
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Figure 1 - continued 



3 RESULTS 

3.1 Dynamics of spatial mixing 

We start the presentation of our results by focussing our at- 
tention on the dynamics of the spatial mixing process. We 
emphasize here that the simulations presented in this paper 
are aimed at understanding the fundamental aspects of the 
dynamics of multiple-population clusters; additional simula- 
tions currently in progress, including a spectrum of masses 
will allow us to more closely compare our theoretical results 
with observations. The panels of Fig[l]show the time evolu- 
tion of the Lagrangian radii of the FG and SG populations 
(along with those of the entire system), and illustrate the 
system's structural evolution and degree of spatial mixing. 
Overlap of all Lagrangian radii for the two populations in- 
dicates complete mixing. As indicated in Fig^ the number 
of initial half-mass relaxation times (of the whole system or 
of the SG subsystem) needed to reach complete mixing of 
the two populations depends strongly on the initial concen- 
tration of the SG subsystem. 

To illustrate the dynamics behind the mixing process 
and the role of relaxation-driven mass loss, Fig. [2] compares 
the time evolution of the ratio of the FG to the SG half- 
mass radii, Rh,fg/ 'Rh,SG, a convenient measure of global 
mixing, for the r25 and the r25isol runs. Initially, the ra- 
tio decreases relatively rapidly, as two-body relaxation acts 
to erase the initial spatial differences between the two pop- 



ulations. However, its rate of change slows as the system 
evolves. The reason for this is that, after reaching core col- 
lapse (at t ~ 17t r h,SG(0) for both the r25 and the r25isol 
systems), the system enters its post-core collapse expansion 
phase. If the cluster did not lose a significant amount of 
mass during the expansion, its relaxation time would in- 
crease, with the half-mass radius and half-mass relaxation 
time growing with time as Rh ~ t 2 ^ 3 and t r h ~ t (see e.g. 
Spitzer 1987). We note that the same behavior, with the 
same implications for the mixing process, is expected no 
matter what mechanism (e.g. three-body binaries, primor- 
dial binaries, stellar evolution mass loss) drives the cluster 
expansion. (See e.g. Gieles et al. 2011 for a recent discussion 
of the role of mass loss due to stellar evolution in provid- 
ing the energy needed to balance the energy flow from the 
cluster half-mass radius.) 

The increasing relaxation time means that the evolution 
of any cluster property driven by two-body relaxation slows 
down accordingly. If we define the cluster dynamical age as 
r — J Q dt/t r h(t) and assume t r h ~ t, it follows that, during 
this phase, r ~ log t and the rate of the cluster dynamical 
aging slows down as dr/dt ~ 1/t. The decreasing rate of 
spatial mixing of the FG and SG populations evident in 
Fig.|2]is a manifestation of the decreasing cluster aging rate. 

In reality, clusters lose mass, and the time evolution of 
trh and r deviate from the expressions derived on the as- 
sumption of zero mass loss. While the spatial mixing rate 
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Figure 2. Time evolution of the ratio of the FG to the SG half- 
mass radii for the r25 (black lower line) and the r25isol (cyan 
upper line) systems. Time is normalized to the initial half-mass 
relaxation time of the whole system in the lower axis, and to the 
initial half-mass relaxation time of the SG subsystem in the upper 
axis. 



decreases with time for both the r25 and the r25isol sys- 
tems, Fig. [5] shows that the evolution of Rh,fg / Rh,SG for 
the r25isol system is significantly slower than for the r25 sys- 
tem. This difference is due to the different mass loss rates 
of the two systems. As the amount of mass lost increases, 
the growth of the cluster relaxation time gradually slows, 
and eventually it starts to decrease, both as a result of the 
decreasing remaining mass and of the eventual contraction 
of the half-mass radius. As the cluster mass decreases so 
does its Jacobi radius, and the outer, less relaxed (and less 
mixed), layers are gradually stripped away. Thus, by enhanc- 
ing the loss of less-mixed layers and slowing the growth of 
the cluster relaxation time, mass loss accelerates the cluster 
dynamical aging process and the evolution toward complete 
spatial mixing. 

We further illustrate these points in Fig. [3] Here 
we focus on the local relaxation time t re iax(r) = 
0.34cr(r) 3 /(G 2 mp(r)\n(0.11N)), which provides a more ac- 
curate measure of the level of mixing expected at different 
distances from the cluster center. Here, a(r) is the cluster 
1-D velocity dispersion, p(r) the cluster mass density and N 
is the total number of stars. Fig. [3] shows the time evolution 
of the local cluster dynamical age r(r, t) — J Q dt/t re i ax (r,t) 
at different distances from the cluster center. The slowdown 
in the time evolution of r when the cluster expansion starts 
(marked by the vertical dotted line) and the differences be- 
tween the dynamical aging of the r25 and the r25isol systems 
are evident in this figure. 

Fig. [3] also shows the radial dependence of the cluster 
dynamical age: as expected, r(r, t) decreases at larger dis- 
tances from the cluster center, and the cluster outer regions 
mix on longer time scales than the inner regions. The im- 
plications of this radial variation of the dynamical age on 
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Figure 3. Time evolution of the cluster dynamical age r (see text 
for definition) for the r25 (solid lines) and the r25isol (dashed 
lines) systems. The five lines for each system refer (from top to 
bottom) to the values of r measured at the 25%, 40%, 50%, 60%, 
75% lagrangian radii. The vertical dotted line marks the time of 
core collapse (see also the top left panel of Fig. 

the mixing process, and its imprint on the radial variation 
of the fraction of SG stars, are discussed in more detail in 
Section [Ql 

The general evolution toward spatial mixing is there- 
fore driven by internal two-body relaxation and accelerated 
by mass loss. The mass loss rate due to two-body relax- 
ation is determined primarily by the strength of the external 
tidal field (see e.g. Vesperini & Heggie 1997, Baumgardt & 
Makino 2003, Gieles & Baumgardt 2008), while the local in- 
ternal relaxation rate is determined by the cluster structural 
properties and their radial variation within the cluster. The 
different processes (and time scales) involved in the mixing 
dynamics imply that, in general, it is not possible to deter- 
mine a universal mixing time scale simply in terms of the 
cluster initial half-mass relaxation time scale. The r25 and 
r25isol systems share the same initial structure but, after a 
given number of initial half-mass relaxation time scales, be- 
cause of their different mass loss history, they have reached 
a significantly different degree of mixing. 

3.2 Dependence on the SG initial concentration 

We have carried out a set of four simulations with different 
SG initial spatial concentrations. In comparing the evolution 
of these systems, we assume they represent clusters with the 
same initial masses and tidal radii (equal to the Jacobi ra- 
dius), but different internal structures determined by the 
concentration of the SG population. As shown in numerous 
investigations based on Fokker-Planck and N-body simula- 
tions (see e.g. Chernoff & Weinberg 1990, Vesperini & Heg- 
gie 1997, Baumgardt &; Makino 2003, Gieles & Baumgardt 
2008), systems with the same tidal radius and mass share 
a similar relaxation-driven dissolution time scale, with only 
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Figure 4. Evolution of Rh,FG / Rh,SG as a function of time nor- 
malized to the half-mass lifetime of each cluster, 7\ / 2 (top panel), 
the ratio of the cluster mass to the initial mass (middle panel), and 
time normalized to a common reference time Tad — 0.17V 0-65 /o; 
(lower panel; see Sect ion [3^2] for a discussion of Tad)- The lines in 
each panel refer to the r25 (black), rlO (green), r5 (purple), and 
the r2p5 (orange) systems. 
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Figure 5. Evolution of Rh,fg/ Rh,SG as a function of time nor- 
malized to T1/2 f° r the r ^ system with (from top to bottom) 
N = 60,000 (top red line), N = 30,000 (middle cyan line), and 
N = 10, 000 (lower black line). 



a weak dependence on the internal structure of the cluster 
(see e.g. Gieles & Baumgardt 2008). The four systems in- 
vestigated here span almost an order of magnitude in initial 
half-mass relaxation time, but only a factor of two in the 
half-mass lifetime (denned as the time needed to lose half of 
their initial mass due to the effects of two-body relaxation) . 
The upper and middle panels of Fig. [4] show the evolution 
of Rh,FG/Rh,SG as a function of time normalized to each 
cluster's half-mass time scale, Ti/2, and as a function of the 
cluster total mass (normalized to its initial value). 

The lower panel of Fig. [4] shows Rh,fg / 'Rh,SG as a 
function of time normalized to a common reference time, 
Tud, denned as T tid = 0.1N°- 65 /uj, where cj = ^GM/3R 3 j. 
The definition of Tud follows from studies showing that this 
functional form provides a good fit to the scaling of the 
relaxation-driven tidal dissolution (or half-mass) time scale 
with initial cluster properties. The exact value of the numer- 
ical factor needed for Tud to match half-mass time scale de- 
pends on the cluster mass spectrum and the cluster structure 
(see e.g. Gieles & Baumgardt 2008 for a discussion on the 
dependence on the latter). Note that Tud is closely related 
to the dynamical family parameter introduced by Chernoff 
& Weinberg 1990 (see their Table 3). 

As just discussed, both internal two-body relaxation 
and (indirectly) mass loss drive the evolution of the FG- 
SG mixing. Despite the initial differences in the SG con- 
centration, all systems approach a state of complete mixing 
(Rh,fg / Rh,SG ~ 1) after approximately the same number 
of Ti/2 and T t %d time scales, and after losing 60-70 percent 
of their mass. 

Fig. [5] compares the evolution of Rk,fg / Rh,SG as a func- 
tion of t/Ti/2 for clusters with the structure of the rlO 

© 0000 RAS, MNRAS 000, 000-000 



Multiple population dynamics 7 




0.5 1.0 2.0 5.0 10.0 20.0 50.0 

t/trh(t) 



Figure 6. Evolution of Rh,fg/ Rh,SG as a function of time t 
relative to the instantaneous cluster half-mass relaxation time 
scale t r h(t), calculated using the 3D half-mass radius, R^. The 
lines refer to the r25 (black), rlO (green), r5 (purple), and the 
r2p5 (orange) systems. 



system but different values of N. The time evolution of 
Rh,FG / Rk,sg shows only a very weak dependence on N. 



3.3 Connection between SG/FG mixing and 

t/Uh(t) 

While, as pointed out at the beginning of this section, the 
goal of the simulations presented in this paper is to under- 
stand the fundamental aspects of multiple-population clus- 
ter dynamics and not to directly compare the results of sim- 
ulations with observations, we nonetheless think it is impor- 
tant to identify those parameters that can be most easily 
measured for real clusters and might be correlated with the 
degree of SG/FG mixing. 

One observable quantity that can be reliably estimated 
for all clusters is the current half- mass relaxation time, 
t r h(t). In Fig. [6] we show the evolution of the degree of mix- 
ing (as measured by Rh,fg / 'Rh,sg) versus the ratio t/t r h(t). 
This figure suggests that t/t r h(t) may be a good observa- 
tional indicator of the extent to which a cluster still retains 
memory of the initial SG segregation. 

It is important to recognize from Fig. [6] that, when in- 
terpreting cluster-to-cluster differences in the current SG- 
FG mixing state, one must also consider possible differ- 
ences in the initial degree of SG concentration. For ex- 
ample, two clusters starting with initial conditions similar 
to our r25 and r5 simulations might be characterized by 
different relative SG-FG spatial distributions even if ob- 
served at the same dynamical phase and after having lost the 
same amount of mass (see Fig. Observations of clusters 
thought to be at the same dynamical phase may thus shed 



light on the extent of the differences in the initial structural 
properties and concentration of the SG subsystem. 

3.4 Spatial mixing and the radial variation of the 
FG-SG number ratio 

As shown in Figs. [I] the FG-SG mixing process occurs more 
efficiently and on a shorter time scale in the cluster inner 
regions, where the local two-body relaxation time scale is 
shorter. In this section we explore further the evolution of 
the radial dependence of the mixing process. 

Fig. [7] shows the radial profiles of the number ratio of 
SG to FG stars Nsg/Nfg for the rlO system at different 
stages of its evolution (the evolution of the radial profile 
shown in this figure is representative of the other systems in- 
vestigated). The system starts with the SG population con- 
centrated in the innermost regions; as the cluster evolves 
and the FG and SG stars mix, the inner flat portion of the 
Nsg/Nfg profiles progressively extends toward the outer 
regions. Complete mixing corresponds to a flat profile ex- 
tending over the entire cluster. 

During the evolution the Nsg/Nfg radial profile is 
characterized by an approximately flat inner region followed 
by a declining outer portion corresponding to regions in- 
creasingly dominated by FG stars. In the outer regions, the 
Nsg /Nfg profile is characterized again by an approximately 
flat portion followed, in some cases, by a weak final rise in 
the profile in the cluster outermost regions. 

As our simulations show, during a large fraction of a 
cluster evolution, the SG and the FG populations are not 
completely mixed and Nsg/Nfg varies with the distance 
from the cluster center. This fact has several implications: 
First, some memory of the initial SG segregation predicted 
by D'Ercole et al. (2008) should still be preserved and ob- 
servable in many clusters today. Second, the observational 
determination of the SG-to-FG number ratio at a given dis- 
tance from the center of a cluster will, in general, differ from 
the global value of that quantity [hereafter we will refer to 
the global value of Nsg/Nfg as (Nsg /N F G) g iob\- Third, 
when exploring cluster-to-cluster differences in the SG-to- 
FG number ratio, we must take into account the possibility 
that part of these differerences might arise if observations 
cover different radial zones of different clusters, relative to 
the half-mass radius. 

The global SG-to-FG number ratio and its cluster-to- 
cluster variation are key ingredients of, and constraints on, 
our ability to understand and model the formation and 
evolution of multiple populations in globular clusters. It is 
therefore important to properly understand the relationship 
between the observational estimates of the SG-to-FG num- 
ber ratio and its actual global value. 

Fig. [8] shows the time evolution of the radius, R g iob, at 
which the local value of Nsg/Nfg equals the global value 
(Nsg /Nfg) glob- The evolution of R g i b is shown only un- 
til t/Ti/2 ~ 1 since, beyond this time, as the system ap- 
proaches complete mixing R g i Q b becomes very noisy and not 
well defined. Fig. [8] shows that Nsg/Nfg generally equals 
(Nsg /Nfg) glob at radii ~ (l — 2)Rh. These results provide a 
first approximate indication of the range of radial distances 
to be targeted in observational studies in order to estimate 
the global SG-to-FG number ratio from observations cover- 
ing only a limited radial range. Additional simulations with 
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Figure 7. Evolution of the radial profile of N$g/Nfg f° r the 
rlO system. Radius is normalized to the instantaneous half-mass 
radius, and N&G /Nfg is normalized to the global SG-to-FG num- 
ber ratio. In each panel the profiles are shown at t/t r h(0) [and 
t/T 1/2 , t/t rh (t)\ approximately equal to 0.2 (0.003,0.2) (black), 
1 (0.015, 1) (blue), 5 (0.08, 3.5) (red), 15 (0.23, 5) (green line), and 
45 (0.7, 10) (magenta). The top panel shows the profiles with a 
log-log scale to highlight the inner structure while the lower panel 
shows the same profiles in linear scale to highlight the outer re- 
gions. 

a spectrum of masses are currently underway, and will al- 
low a refined estimate of R g iob and more direct comparisons 
with observations. 



3.5 Evolution of the global SG-to-FG number 
ratio 

As shown by D'Ercole et al. (2008), the transition from a 
cluster initially dominated by FG stars to one with sim- 
ilar numbers of SG and FG stars (or even dominated by 
SG stars) occurs mainly during the cluster's early evolu- 



in 
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Figure 8. Time evolution of the ratio of the radius R g i bi 
at which the local value of N$g/Nfg equals the global value 
{NsG/NFG)globi to the cluster half-mass radius, R^. Black, 
green, purple and orange lines refer, respectively, to the r£5, rlO, 
r5, and r2p5 simulations. 

tion. This transition is due to the loss of FG stars during 
the cluster's early expansion triggered by the loss of SNII 
eject a. 

During the subsequent long-term evolution, mass loss 
due to two-body relaxation removes both FG and SG stars 
from the cluster. Until the two populations are completely 
mixed, two-body relaxation will still lead to a slight pref- 
erential loss of FG stars and further increase the SG-to-FG 
number ratio (although this increase is much smaller than 
that during the cluster early evolution; see D'Ercole et al. 
2008). Fig. [9] shows the time evolution of the global SG- 
to-FG number ratio, (Nsg /NFG) g iob- For the systems stud- 
ied here, during the long-term evolution (Nsg /NFG) g iob in- 
creases by only a factor of ~ 1.3 — 1.4. 

If stars with chemical properties typical of observed SG 
populations can form only in globular clusters, the SG stars 
observed in the Galactic halo (see e.g. Martell & Grebel 
2010, Carretta et al. 2010b) must have been lost by clusters 
mainly during the long-term evolution phase explored in this 
paper. As discussed by Vesperini et al. (2010; see also Car- 
retta et al. 2010b, Schaerer & Charbonnel 2011), the fraction 
of SG stars in the Galactic halo is a very important quan- 
tity that can shed light on the connection between globular 
clusters, their formation and dynamical history, and their 
contribution to the assembly of the Galactic halo. 



4 DISCUSSION AND CONCLUSIONS 

Our previous study of the formation and evolution of multi- 
ple populations (D'Ercole et al. 2008) showed that SG stars 
form in the innermost regions of a cluster initially dominated 
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Figure 9. Time evolution of {Nsg/Nfg) glob f° r the simulations 
r25 (black), r7#(green), r5(purple), r2p5 (orange). 



by FG stars. Following the early cluster expansion triggered 
by mass loss due to SNII, a large fraction of FG stars escape 
the cluster, leaving a system with comparable numbers of 
FG and SG stars, with the SG stars strongly segregated near 
the cluster center. The simulations presented here follow the 
subsequent long-term evolution of these clusters, driven by 
two-body relaxation. By means of a survey of N-body sim- 
ulations, we have investigated the dynamical evolution of 
multiple-population clusters, focusing our attention on clus- 
ter structural evolution and the FG-SG mixing process. We 
have studied initial conditions characterized by different de- 
grees of the initial concentration of the SG subsystem. 

Two-body relaxation is the key process driving the mix- 
ing of the two populations. As the local two-body relaxation 
time scale increases in the cluster outer regions, this process 
becomes less efficient at larger distances from the cluster 
center. Moreover, as soon as an energy source in the cluster 
core causes a cluster to start expanding, the cluster relax- 
ation time increases with time and the rate of its dynamical 
aging increasingly slows down. A significant slowdown in the 
evolution toward complete mixing ensues (see Figs. [2] and [3]). 

Our simulations show that as cluster evolution contin- 
ues, relaxation-driven mass loss affects the mixing process 
by slowing and eventually reverting the growth of the clus- 
ter relaxation time, and by causing the loss of the cluster's 
outer unmixed layers. By comparing the evolution of clus- 
ters with the same initial structure but different mass-loss 
rates, we have quantified this process and illustrated the 
joint role of internal relaxation and tidal mass loss in the 
evolution toward complete mixing. Our simulations imply 
that, unless a cluster has lost a significant fraction of its 
mass ( 60 — 70%; as discussed in the paper, we refer here 
to the mass lost during the cluster long-term evolution, and 
not to the early loss of FG stars), FG-SG mixing will not 



be complete and the SG stars will be centrally concentrated 
relative to the FG population (see Fig.|4|. 

To further explore the radial dependence of the mix- 
ing process, we have followed the time evolution of the ra- 
dial profile of the SG-to-FG number ratio, Nsg/Nfg (see 
Fig.[7|. We have shown that as SG-FG mixing proceeds, the 
Nsg/Nfg profile can be divided into three different regions: 
(1) an inner region where Nsg/Nfg is fiat; (2) an interme- 
diate region where the FG is increasingly dominant, with 
Nsg/Nfg declining with increasing radius; (3) an outer re- 
gion where the Nsg/Nfg profile flattens (the outermost 
cluster regions are in some cases characterized by a slightly 
rising Nsg/Nfg profile). 

Since, as discussed above, reaching complete mixing 
(and therefore a flat Nsg/Nfg radial profile) requires the 
cluster to be in an advanced phase of dynamical evolution, 
with the loss of a significant fraction of stars, we expect that 
a Nsg/Nfg profile characterized by the regions described 
above will be found in all clusters which have not undergone 
strong relaxation-driven mass loss — possibly a large fraction 
of the Galactic globular cluster population. Although a de- 
tailed and quantitative comparison with observations is be- 
yond the scope of the analysis carried out in this paper, 
our study provides a preliminary indication of the link be- 
tween the level of SG/FG mixing and parameters that can 
be determined observationally for real clusters. Specifically, 
we have illustrated the relation between SG/FG mixing and 
the ratio t/t r h(t) (see Figj6|. 

A few observational studies have explored the differ- 
ences in the spatial distribution of FG and SG stars in glob- 
ular clusters, and find that, in agreement with the prediction 
of the formation and evolution models presented in D'Ercole 
et al. (2008) and the results presented in this paper, SG 
stars do tend to be concentrated in the cluster inner regions 
(Bellini et al. 2009, Carretta et al. 2010a, Lardo et al. 2011, 
Kravtsov et al. 2010, 2011 Nataf et al. 2011, Johnson & Pi- 
lachowski 2012, Milone et al. 2012). It is interesting to note 
the similarity between the overall shape of the Nsg/Nfg 
profile found in our simulations and the observational pro- 
files as presented in Lardo et al. (2011), Bellini et al. (2009), 
and Milone et al. (2012). It is interesting to point out that 
all these clusters have values of t/t r h(t) (assuming a com- 
mon age of t — 11.5 Gyr and t r h calculated from the Harris 
(1996, 2010 edition) catalogue, taking into account the ap- 
proximate relation between the 2D half-mass radius Rh,2D, 
and the 3D half-mass radius Rh, Rh — |^,2d) such that 
the presence of a radial gradient in the Nsg/Nfg profile is 
expected, and in agreement with the simulations presented 
in this paper (see Fig|6]). 

Unless a cluster is completely mixed, the SG-to-FG 
number ratio determined from observations at a given dis- 
tance from the cluster center is in general different from the 
global value. This implies that care must be used in com- 
paring the SG-to-FG number ratio in different clusters: if 
different clusters are observed at different distances from the 
center, cluster-to-cluster differences in this number might be 
found even for systems with the same (Nsg /NFG) g iob- Our 
simulations indicate that, so long as mixing is not complete, 
the local value of Nsg/Nfg measured at R m (1 — 2)Rh is 
approximately equal to the global SG-to-FG number ratio. 
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